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Strange Magnetic Moment of The Nucleon from Lattice QCD 
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We calculate the strange magnetic moment of the nucleon on a quenched 16^ x 24 lattice aX (3 — 6.0, and with 
Wilson fermions at k = 0.148, 0.152, and 0.154. The strange quark contribution from the disconnected insertion 
is estimated stochastically by employing the Z2 noise method. Using an unbiased subtraction along with the help 
of charge conjugation and hermiticity, we reduce the error by a factor of 2 with negligible overhead. Our result 
is Gli = -0.28±0.10Aiiv. 



1. INTRODUCTION 

Recently there has been a considerable inter- 
est on the strange magnetic moment of the nu- 
cleon (G^(0)) as it can provide important in- 
formation about the internal quark structure of 
the nucleon There are numerous theoretical 
predictions over a wide range in between —0.75 
to -1-0.30/^ AT. Most of them are negative and are 
in between —0.45 to —0.25 (for summaries see 
j2||^). On the other hand, recent experiments 
of SAMPLE collaboration on parity violating 
electron scattering [|[|| suggest that the strange 
magnetic moment of the nucleon could be posi- 
tive. Two of their recent results are G%j(Q^ = 
O.lGeV^) ^ +0.23 ± 0.37 ± 0.15 zbO.ig^iAr g and 
+0.61 ± 0.17 ± 0.21 ± 0.l9fiN !]• However, com- 
bined with data from the deuteron, the new result 
is closer to zero with a relatively larger error 

Previously, our group studied the strange quark 
magnetic moment on a quenched lattice and ob- 
tained G%i{0) = -0.36± 0.20^*^ II . We extend 
that calculation and improve the results by em- 
ploying CH-transformation and unbiased subtrac- 
tion method, to be discussed subsequently. Fi- 
nally, we will present our new results on G%j{0). 

2. FORMULATION 

We follow the same lattice formulation of the 
electromagnetic vector current as in Rcf. MM- 



The magnetic moment of the nucleon is extracted 
from the ratio of the three and two-point correla- 
tion function. Two point function is defined as 

G%%{t,p) = 

Ve-^^->|T(x"(x)x''(0)) |0), (1) 



and the three point function is 
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(0|T(x"(x/)J4a:)x'^(0)) |0). (2) 

Here, is the point-spilt electromagnetic vector 
current as defined in Ref . Q . Finally, taking the 
ratio 
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where G\,i{q^) is the magnetic moment of the 
nucleon on lattice at four-momentum transfer 
square . 

The three point function incorporates both 
connected insertion (CI) and disconnected inser- 
tion (DI). CI arises due to valence and connected 
sea quarks and DI originates from the vacuum 
polarization of the disconnected sea quarks 0. 
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For the nucleoli magnetic moment, strange quarks 
current, 57^ s, contributes only through the DI 
(Fig. 1), and so, we will concentrate only on DI 
contribution of Eq. (3) . 
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Figure 1. Disconnected insertion diagrams to 
show the strange quark contribution to the nu- 
cleon magnetic moment. 



3. LATTICE CALCULATION 

In general, DI is very difficult to calculate as 
it involves loop effect from the self-contraction of 
the sea quarks of the current. As in the previ- 
ous DI calculations [p|-p^, for Eq. (3), the cur- 
rent insertion time (t) from the nucleon source 
to its sink is summed over, and also the aver- 
age over three spatial directions is considered. 
With this average, this ratio R(3pt,2pt)(Eq. (3)) 
becomes proportional to constant -I- tf[\q\/Eq + 
™]^a/(9^) 0- Therefore, by extracting the 
slope of R(3pt,2pt) vs tf plot, one can evalu- 
ate G^{q^). As shown in Fig. 1, the three 
point function of DI consists of two parts : nu- 
cleon propagator (two-point function) and a loop 
contribution. Similar to other DI calculations 
|8|-[l0|], to evaluate the trace for the quark loop, 
we use Z2 noise trace estimator with a stochas- 
tic algorithm. It has been shown earlier that by 
employing CH-transformation and unbiased sub- 
traction method one can reduce the error in DI 
substantially ||9|,|l0||. Similar error reduction pro- 
cedure is adopted here. By CH-transformation 
we mean the charge conjugation (C) and Eu- 
clidean hermiticity (H). Under charge conju- 
gation, the fcrmion matrix (M) transforms as 
M-^{x, y, U) = C-^M-^{y, x, U*)C. Along with 



this, using hermiticity property of the fermion 
matrix {j5M~^{x, x + a^)^^ — M^^^x + a^, x)), 
one can show that under CH-transformation, 

M-\x,y,U) = C-^^^M-^* {x,v,U*)-i5C. (4) 

Using this CH-transformation, it is possible to 
show that the three point function is either real 
or imaginary, and for the vector current it is 
real. Using only Euclidean hermiticity, one can 
also show that the loop part is imaginary, and 
thus, after Fourier transformation it becomes 
real. Therefore, to get the three point function 
one needs to multiply only real part (real with 
Fourier transformation) of loop with the real part 
(real with Fourier transformation) of the two- 
point function. Thus, one can avoid unneces- 
sary imaginary x imaginary multiplication, and 
thereby can eliminate unwanted noises. 

Along with this CH-transformation, as in angu- 
lar momentum calculation | ]To| , we also adopt an 
unbiased subtraction method (for details discus- 
sion about unbiased subtraction, see Ref. [p|,p^). 
The trace of an inverse matrix can be esti- 
mated as : 



E 



(5) 



where 77s are Z2 noise vectors, O^s are a set of 
traceless matrices, and ApS are variational coef- 
ficients which can be tuned to get the minimum 
variance in trace estimation. It was found earlier 
P that a set of traceless matrices obtained from 
the hopping parameter expansion of the propa- 
gaptor 

= / 



kD 



can offset the off-diagonal contribution to the 
variance of the trace, and thus they can be used 
as a possible choice of O*-^' matrices. For our 
calculation we use traceless matrices up to order 
D\ 

Numerical calculation is done on a quenched 
16'^ X 24 lattice at /? = 6.0 with Wilson fermions, 
and for 100 configurations. Results are obtained 
for relatively light quarks with k = 0.148,0.152, 
and 0.154. For the stochastic trace estimation, 
300 noises per configuration is used. 
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Figure 2. The summed ratios (R(3pt,2pt)) of Eq. 
(3) are plotted for different final time-slices. Re- 
sults with and without subtractions (along with 
CH-transformation) are represented by solid and 
open circles, respectively. Ratios with subtrac- 
tions are shifted slightly towards right. 



Figure 3. Monopole (top) and dipole (bottom) 
fitting for the strangeness magnetic form fac- 
tor G%[{q^) of proton. Fitted G%i{Q) values are 
shown by o. 



4. RESULTS 

As mentioned earlier, the strangeness content 
of the nucleon magnetic moment , G|^(0), ap- 
pears only in the DI. The numerical procedure 
for extracting G\j{q^) has been given in detail in 
Ref. 0. First, we fix the sea quark mass (loop 
mass) at the strange quark mass (at k — 0.154), 
and then, for various valence quark masses (k = 
0.148,0.152, and 0.154), we calculate slopes for 
the summed ratios, R(3pt,2pt), with respect to 
various final time slices. In Fig. (2), we plot 
these summed ratios as a function of time slices 
(at \q\ = 27r/L). As mentioned earlier, the slope 
of this plot provides the DI contribution of the 
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magnetic moment, G\j{q^). One can easily no- 
tice from these plots that we reduce the error-bar 
of previous calculation (ratio with open circle) 
substantially (by a factor of ~ 2). These ratios 
are calculated at various g^, and finally, to ob- 
tain the physical G%j(q^), we extrapolate the va- 
lence quark mass to the chiral limit with a form 
A + B^fh + TTis, where fh is the average u and 
d quark mass and rus is the strange quark mass, 
respectively. This form of mg-dependence is cho- 
sen as G%j is proportional to the kaon mass, m^. 
Finally, to obtain the strange quark contribution 
to the magnetic moment, G^j(O), we fit different 
G^(g^) values with monopole and dipole forms, 
as shown in Fig. (3). For all these fittings (Fig. 
2 and 3), required correlations among different 
quantities are taken into account by covariant ma- 
trices, and all errors are estimated from jackknif- 
ing the fitting procedure. 

From the monopole fit G\,j{q^), we extract 
G^(0) = -0.28 ± O.lO^iAr, with the monopole 
mass, niM/mN = 0.80 ± 0.09. It should be re- 
called that the previous result (before subtrac- 
tion) for G%j{0) was -0.36 ± 0.20^*^ §, i.e., 
we have reduced error by almost a factor of 2. 
Another fit of G%j{q'^) with dipole form yields 
G^,j(0) = —0.24 ±0.08^ AT, and the corresponding 
dipole mass is rriD/mM = 1.44±0.16. As argued 
in Ref. 0J|,|l^, we choose the monopole form 
over dipole for DI. Some other useful quantities, 
e.g., the total sea quark contribution, G|f''(0), are 
also calculated. Taking account CI results from 
Ref. [Q, we also evaluate the proton and neutron 
magnetic moments. Since the total sea quark con- 
tribution, G||°(0), is small, these moments are al- 
most same as in Ref. . In Table 1 , we summarize 
our result for various quantities. A comparison is 
also made to show the error reduction. 

5. CONCLUSIONS 

We reduce the error of strange quark magnetic 
moment, as estimated previously Q, by a fac- 
tor of 2. Strangeness magnetic moment of the 
nucleon is still negative and is —0.28 ± 0.10/xjv. 
Once more, we prove that the error reduction 
procedure, by implementing CH-transformation 
and unbiased subtraction method, works well in 



Table 1 



Nucleon magnetic moment (in /ijv) before and af- 
ter error reduction. 





Before 


After 


G (0)77107^0 


-0.36(20) 


-0.28(10) 


G%[ifi)dipo 


-0.27(12) 


-0.24(8) 


M \^ J mono 


-0.65(30) 


-0.60(16) 


Gm iS^^dipo 


-0.45(16) 


-0.40(12) 


Gf/(0) 


-0.097(37) 


-0.106(28) 


GM{fi)proton — 


2.53(10) 


2.52(8) 


GM{0)neutron — 


-1.71(12) 


-1.71(9) 


GAf(0)p,o ^, 

GAf(0)„.„ - 


-0.68(7) 


-0.68(5) 



reducing error in the disconnected insertion. Fi- 
nally, we like to mention that the finding of this 
work is in disagreement with the results of the 
next talk [|l|. 
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